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Abstract. We present preliminary results of terrestrial planet formation using on the one hand 
classical numerical integration of hundreds of small bodies on CPUs and on the other hand—for 
comparison reasons—the results of our GPU code with thousands of small bodies which then 
merge to larger ones. To be able to determine the outcome of collision events we use our smooth 
particle hydrodynamics (SPH) code which tracks how water is lost during such events. 
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1. Introduction and scenarios 

The outcome of numerical simulations of planet formation is extremely sensitive to 
initial conditions like disk surface density, total mass, the initial distribution of plan¬ 
etesimals and planetary embryos after the gas phase in the protoplanetary disk, and the 
dynamical model. Hence, the resulting planetary systems and even statistical statements 
are highly biased by the choice of these parameters. While many results on formation of 
planets and their wat er content have been achieved, several open questions remain (e.g., 
Raymond et al. 2014 1. 

Here, we present preliminary dynamical results for a model including two gas giants, 
Jupiter and Saturn. By SPH simulation of collision events we can deduce the outcome of 
collisions in terms of merging/fragmentation and water loss d epending on mass and siz e 
of the bodies, collision velocity, and impact angle (Fig.[5jD and iMaindl et al.ll2013ll2M^ . 
which addresses the yet unclear co ntribution of includ ing water transport probabilities 
in planet formation simulations Icf. Ilzidoro et al.ll2014ll . __ 

We take the outcome of the last phases of the Grand Tack scenario ( Walsh et al. 201 ih 
as initial conditions for our simulations. Jupiter and Saturn are in their actual positions 
in the Solar System and we distribute small bodies in the mass range 3 x 10“® Mq < 
m < 3 X 10“^ M© and semi-major axes 0.4 AU < a < 2.7 AU with small inclinations and 
small eccentricities. Bodies initially outside 2 AU are assumed to have a water content of 
10 percent. 


2. Preliminary simulation results 

Figure [1] shows a typical result out of 30 different runs after an integration time o f 
5Myrs using the Lie integration method with adaptive step size ( Eggl fc DvorakIboiOll . 
We are able to include thousands of gravitating bodies instead of hundreds by using 
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Figure 1. Eccentricities and inclinations for a typical scenario. The water content is color-coded, 
circle sizes denote the bodies’ masses after 5 Myrs. The largest bodies are of ~ 30 Lnnar masses. 




Figure 2. a) The GPU-over-CPU computational gain first grows with the number of bodies 
involved and settles at a factor of approximately 100. b) Low-resolution SPH simulation of two 
Ceres-mass bodies colliding at Ikm/s and a 30° impact angle (50k SPH particles). Blue dots: 
30wt-% water ice mantle on the target, black dots: solid basalt projectile. 


our n ew massively parallelized GPU code (a descendant of the one described in ISiilil 
2013h which is up to two orders of magnitude faster than our CPU n-body integrator (cf. 


Fig. [5^). This large number of bodies will take care of dynamical friction with a precision 
not achieved up to now (e.g., O’Brien et al.ll2006l) . 


3. Conclusions 

While th ere exist several results from differ ent groups with different underlying assump¬ 
tions (e.g., Izidoro et al.l 20l4 HansenI 20091) . we believe that independent computations 
including phenomena like dynamical friction and realistic collisional water transfer are 
desirable for answering the key questions of terrestrial planet formation. 
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